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The partition function of the two-dimensional Ising model is exactly obtained on a 
lattice with a twisted boundary condition. The continuum limit of the model off the 
critical temperature is found to give the mass-deformed Ising conformal field theory 
(CFT) on the torus with the complex structure r. We find that the renormalization 
group (RG) flow of the mass parameter can be holographically described in terms of 
the three-dimensional gravity including a scalar field with a simple nonlinear kinetic 
function and a quadratic potential. 


1. Introduction 

The AdS/CFT correspondence[1, 2, 3] has been giving valuable informations for various field 
theories especially in the strong coupling regions. In particular, the correspondence between 
the classical gravities and the large N (gauge) field theories has been most intensively studied. 
Among many tests on this correspondence, the exactly soluble models in low dimensions, 
say, the two-dimensional conformal field theories, are expected to shed light. The relation 
of AdS gravity with the Virasoro algebra of CFT is based on the pioneering work of Brown 
and Henneaux[4], and the minimal conformal field theories in two dimensions have been 
discussed in relation to the AdS gravity in three dimensions[5, 6]. Subsequently extensive 
analyses have been done for the Wn minimal conformal field theories and many evidences 
are presented for the relevance of the higher spin field theories together with gravity for these 
minimal models in the large N limit[7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. 

As for the recent progress in the stronger statement of the AdS/CFT correspondence, that 
is, the correspondence between the “quantum” gravity or the string theory and the finite N 
(gauge) held theory, the finite N effect of the one-dimensional supersymmetric gauge theory 
with 16 supercharges (the finite N BFSS matrix model[23]) has been directly examined using 
the Monte Carlo simulation[24, 25], where the a'-corrections to the Type IIA supergravity 
[26] are reproduced from the gauge theory. This strongly motivates to regard held theories in 
some category as candidates of the quantum gravity at least in the semi-classical meaning. 
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In fact, in Ref. [27], it has been proposed that the two-dimensional minimal models are 
holographically dual to three-dimensional quantum gravities. In particular, the simplest 
minimal model, namely the Ising conformal field theory, is conjectured to be dual with the 
three-dimensional Euclidean pure quantum gravity. 

We here briefly review the argument in Ref. [27]. Let us consider the three-dimensional 
quantum gravity with negative cosmological constant. The “quantum” here means that we 
integrate over all the possible three-dimensional metric. In this path-integral, the boundary 
of the three-dimensional geometry is fixed to the 2 -torus with the complex structure r, and 
thus the partition function of the gravity is a function of r and f. The important fact is 
that the smooth classical solutions of the three-dimensional pure gravity are restricted to the 
locally AdS geometries, which are obtained by the SL( 2 , Z) transformation from the thermal 
AdS geometry. Thus, if we can use the semi-classical approach, the partition function of the 
quantum gravity is obtained by summing up the classical contributions of each solution with 
some quantum corrections. The expansion parameter of the saddle point approximation is 
the inverse of the quantity, 


c = 


3 L 

2GV 


( 1 . 1 ) 


where L is the AdS radius and Gat is the Newton constant. The quantity c equals to the 
central charge of the two-dimensional Virasoro algebra which appears at the boundary of 
the asymptotically AdS geometry [4], Therefore the senri-classical approximation is usually 
expected to work only for large c region. The most important assumption in Ref. [27] is that 
the partition function can be evaluated by summing up the classical geometries even in the 
strong coupling region, c ~ 1. Once this assumption is accepted, the partition function can 
be written as 


^grav(r,r)= - z 'vac(7T7'r), (1-2) 

7£SL(2,Z)/r c 

where T c is a subgroup of SL{ 2, Z) which does not change the topology of the three- 
dimensional geometry and Z vac (r, r) is the contribution from the thermal AdS geometry. 
The points are that Z vac (r, f) can be calculated explicitly and that the symmetry T c is 
enhanced at specific values of c < 1 where the summation of ( 1 . 2 ) becomes a finite sum. 
Note that these values exactly equal to the central charges allowed for the two-dimensional 
minimal models. In particular, when c = 7 , the partition function of the gravity (1.2) repro¬ 
duces that of the two-dimensional c = 7 conformal field theory, that is, the Ising conformal 
field theory. This result suggests the duality between the c = 4 minimal model and the 
three-dimensional quantum pure gravity. 

This argument actually gives the correspondence between the c = ^ conformal fixed point 
of the theory space of the two-dimensional field theory and the three-dimensional quantum 
gravity. Then it is natural to expect that there will be a gravity description at least in the 
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vicinity of the conformal fixed point in the theory space of the two-dimensional quantum field 
theories. In the context of the AdS/CFT correspondence, the renormalization group (RG) 
flow of a coupling constant in the field theory can be identified with the classical trajectory 
of the corresponding bulk field in the asymptotically AdS geometry, which is called the 
holographic RG [28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38] (see also [39]). In this scheme, 
the radial coordinate of the AdS gravity can be identified with the RG parameter and the 
AdS boundary corresponds to the conformal fixed point. In the case of the Ising model, 
although the conformal fixed point corresponds to the “quantum” gravity, we can expect 
that the central assumption in Ref. [27], namely, only the classical solutions contribute to 
the partition function of the quantum gravity, still works at least in the vicinity of the AdS 
solution even after adding an additional field in the bulk theory. Then we can examine the 
RG structure around the c = 1/2 conformal fixed point using the traditional technique of 
the holographic RG. 

The purpose of our paper is to work out the exact solution of the Ising model partition 
function on the lattice corresponding to the torus with the generic complex structure r, 
and study the corresponding holographic RG structure of the continuum theory around the 
conformal fixed point. Usually the partition function of Ising model in two dimensions is 
obtained for rectangular lattices with the periodic boundary condition, and one obtains the 
Ising conformal field theory on a rectangular torus (r = i) by taking the continuum limit at 
the critical temperature [40, 41, 42]. However, the most general torus has a complex structure 
with the parameter r, representing the shape of the torus. In the literature, we have found 
no explicit solution of the Ising model partition function on the lattice corresponding to 
the torus with the generic complex structure, although there have been works to obtain 
hnitized conformal spectrum of Ising model on manifold of various topology using corner- 
transfer matrix, Yang-Baxter technique, or thermodynamic Bethe Ansatz[43, 44, 45, 46]. We 
explicitly compute the partition function of the Ising model on the twisted lattice, that is, a 
lattice with such a boundary condition that the “space” position is shifted to some amount 
when one goes around the “time” direction, and show that it ends up with the torus with 
the complex structure in the continuum limit. By taking the continuum limit at off-critical 
temperature with an appropriate scaling, we identify the deviation parameter as the mass of 
the Ising conformal field theory on the torus with the complex structure. We also work out 
the classical solution of the three-dimensional Einstein gravity with a single scalar field. We 
find that a simple nonlinear kinetic function and a simple quadratic potential for the scalar 
field can capture the RG flow from the Ising field theory with the central charge c = 1/2 at 
the ultraviolet fixed point towards the c = 0 case at the infrared by using the technique of 
the holographic RG via the Hamilton-Jacobi equation of the gravity. 
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This paper is organized as follows: In the next section, the partition function of the two- 
dimensional Ising model is obtained on the twisted lattice representing the discretized version 
of the torus with the complex structure. In section 3, the continuum limit of the partition 
function is obtained retaining the deviation from critical temperature, which results in a 
mass term for a free Majorana fermion. In section 4, the holographic description of the 
renormalization group flow is worked out for the massive Majorana fermion in terms of the 
Hamilton-Jacobi equation. Section 5 is devoted to a summary of our results and a discussion. 
Some technical details in computing the partition function of the Ising model is summarized 
in Appendix A. Some details of partition function of two-dimensional massive Majorana 
fermion on the torus is given in Appendix B. 

2. Partition function of the 2D Ising model on the twisted lattice 

Let us consider the 2D Ising model on a rectangular lattice with the size n x m. There 
is a ’’spin” s(x) = ±1 at each site x = {x\,X 2 ) (x\ = 1 ,..., n, X 2 = 1, ■ ■ ■, m) and the 
Hamiltonian of the system is given by 

H = - Ji s(a:)s(x + 1 ) — J 2 ^ s(x)s(x + 2 ), ( 2 . 1 ) 

X X 


where J\ (J 2 ) and 1 (2) are the coupling constant and the unit vector in the “space” (“time”) 
direction, respectively. As for the boundary condition, we impose the periodic boundary 



Fig. 1 The spin degrees of freedoms are on the sites of the n x m lattice. In our case, we 
take the usual periodic boundary condition for the x\ direction but we adopt the twisted 
boundary condition for X 2 , namely, we identify ( x\ + p,x 2 + m) with (aq,^)- 


condition to the space direction and the twisted or shifted boundary condition with an 
integer parameter p € Z to the time direction, namely, the space position is shifted by p 
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when one goes around the time direction (see Fig. 1): 

s(xi + n,x 2 ) = s(xi,X 2 ), s(x\ + p,x 2 + m) = s(xi,X 2 ). (2.2) 

In the following, we evaluate the partition function, 

^ J2 e ~ PH > ( 2 - 3 ) 

{s(a:)=±l} 

under this boundary condition. 

Let us first define the matrices with the size of 2 n , 

S{ = 1 2 <8> ■ ■ ■ (8> cr { <g> ■ ■ ■ <g> 1 2 , (2.4) 

where the Pauli matrices a 1 ( i = 1,2, 3) are placed at the k-th position (k = 1,..., n). Using 
these matrices, the transfer matrix of this system is written as 

T = (2sinh2a) ¥ 1414, (2-5) 


with 

n n 

V a = W exp (aSi) , V b = exp ( bS%Sl +1 ) , 
k =1 k =1 

where a and b are given by a = /3 J\ and b = (3 J 2 and a is defined by 

sinh 2 a = 


( 2 . 6 ) 


(2.7) 


sinh 2 a 

Note that the system is in the ordered phase for a < b and is in the disordered phase for 
a > b. We further define the ’’shift matrix” £ whose components are explicitly given by 


( 2 . 8 ) 


'Sk,2k-l ~ ^2 n_1 +fc,2fc = 1 , (k — 1 , • • • , 2 n : ) 

0 , (others) 

We see that £ has the property, 

£-^£ = Si_ l7 (S' n+1 = S{). (2.9) 

We can then write the partition function in Eq.(2.3) with the boundary condition in Eq.(2.2) 


as 


Z = TV (T m £ p ). (2.10) 

In order to estimate the transfer matrix, we here construct the Dirac matrices of Spin(2n ) 


as 


r 2 fc-! = Si52 ■ • • Sfc-iSfc, T 2k = SlSl--S 1 k _ 1 S 2 k , (k = 1,... ,n) (2.11) 

and 

U=SlSl---S'=i n T 1 ---T 2 n, ( 2 . 12 ) 
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which satisfy the Clifford algebra, 


{IV, r„ } = 2S lhV , {r M ,C/}=0, U 2 = 1. (fi,u = !,-■■ ,2 n) (2.13) 


As usual, the generators of Spin{2n ) in the chiral and anti-chiral representations are defined 
as 

j± = J[r M ,r v ]c/ ± , (2.14) 

where 

U± = ^(1 ± U). (2.15) 

Using them, we can also divide V a , 14 and E into the chiral and anti-chiral sectors as 
V± = 14 U±, V± = 14 U± and E± = E U±. In the following, we use the notation, 

H ± = (v±) 1,2 v±(vtf /2 ■ (2.16) 

For E±, we can easily see 

(U a ± )- 1/2 E ± (U a ± ) 1/2 = E ± . (2.17) 

Then the partition function Eq.(2.3) can be written as 

Z =(2sinh2 a)^ TV (H™YF + U+ + H™YF_U-) 

=(2sinh2a)^{Tr + (77™E*{) + TV_ (H m T p _)\, (2.18) 


where TV-|- denote the trace over the chiral and anti-chiral sectors of the spin representation 
of Spin(2n). respectively. 

In evaluating (2.18), the following fact is useful: Suppose that 3 A S Spin(2n\<C) in the 
fundamental representation is transformed into the “canonical form” by T £ 0(2n ) as 

n n 

T t AT = 0 = 0 R(-iO k ), (0 k € C) (2.19) 

k=l k=l 


where J /ijy are the generators of Spin(2n; C) in the fundamental representation and R(9) is 
the two-dimensional rotation matrix with the (complex) angle 9, 


R(0) = 


( cos (9) 
ysin(0) 


— sin(0) 
cos(0) 


Then we can write Tr±(A) as 


Tr±(A) 


1 

2 



n 

± det(T) 2sinh 
fe=l 



( 2 . 20 ) 


( 2 . 21 ) 


Thus, we first express 77™E({ in the fundamental representation and then transform them 
into the canonical form using appropriate matrices T± € 0(2n). 
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The computation along this strategy is straightforward and we summarize it in Appendix 
A. As a preparation to show the result, we define the quantity 77 (/ = !,■•• , 2 n) as the 


positive solution of the equation, 



cosh 7 / = cosh 2 a cosh 26 — cos f— ) sinh 2 d sinh 26, 

V n ) 

( 2 . 22 ) 

and 

7T pi 

7/ - 7/ 2 

(2.23) 


mn 


Note that 7 j satisfy the relation: 

e m72_ e m7i*_ 

(2.24) 


In addition, we will often use 70 = 72 n in the following. Combining the above consideration 
and the results (A24) and (A25), we obtain the partition function of the two-dimensional 
Ising model in the ordered phase with the twisted boundary condition in Eq.(2.2); 

\ mn 4 

Z = — (2 sinh 2a) ~ ^ ' Z,;, (2.25) 

i —1 

with 


r= 1 


r—1 


r—1 


Zi = Ri n 12 cosh (y 72 r- 

r=l 

[f] 

Z 2 = R 2 12 sinh {-^l 2 r 

P 3 2 cosh 


'J e 2 72fc_1 ) Pi, 
Vfc=i / 

^ n \ 

ne>- p 2 , 

<fc=l / 


(2.26) 


/m 

Z 3 = 2 cosh —70 


V 2 


m 


r=l 
r n —11 


pm 

[-2 72r 


] ■ (! + e _m7 °) P 3 , 


Oo=l 


Z 4 = 2 sinh ( —70 


P 4 12 sinh 


r =1 


m, 


-72 r 


J e? 7 - • (1 - e- m7 °) P 4 , 


^fc=i 


where we have used the reflection property (2.24) and introduced 


Pi = < 


P 3 = < 


1 (n : even) 

2 cosh (?r7n) ( n : odd, p : even) , P 2 = < 

2 sinh (y7n) (p '■ odd, p : odd) 

2 cosh (y7n) (n ■ even, p : even) 

2 sinh (y7n) (ra : even, p : odd) , P 4 = < 

1 (n : odd) 


1 (n : even) 

2 sinh (y 7 n) (w : odd, p : even) 

2 cosh (y 7 n) (^ : odd, p : odd) 

2 sinh (y 7 n) (n, : even, p : even) 

2 cosh (y 7 n) (^ : even, p : odd) 

1 (n : odd) 

(2.27) 
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and 


Pi = 

(nil + e-™^- 1 ! 2 ) 

r =1 


P2 = 

[?] 

r—1 

(l - (-l) P 7- 1) ».- 1 e-"* T ") 

P 3 = 

[V] 

( n |l + e _m ^| 2 ) 

r= 1 


Pi = 

[V] 

(n |i ~ e - m ^i 2 ) 

r= 1 

(l-(—1)7—”*>). 


The result Eq.(2.26) constitutes our new result of partition function of the two-dimensional 
Ising model on a twisted lattice that gives the discretized version of torus with the complex 
structure t 

T) Tfl 

T = T\ + zr 2 , T\ = — , r 2 = — . (2.29) 

n n 

whose continuum limit gives the torus with the complex structure t as described in the next 
section. 


3. Continuum limit of the partition function 

We next consider the continuum limit, that is, the limit of m, n,p —>• oo with fixing the ratios 
ri, t 2 in Eq.(2.29). In taking this limit, we also tune the coupling constant so that the theory 
properly reaches to a continuum theory around the critical point. For simplicity, we consider 
the case of J\ = ,/ 2 , that is, a = b = K. Under this condition, we defined the parameter [i 
through the relation 

= sinh2 K +---^ - 1. (3.1) 

4n 2 2 \ sinh 2K J 

Recalling that the critical temperature K* is given by sinh2iv* = 1, the parameter /i 
expresses a deviation from the critical temperature for a finite n [42]. Note that we take the 
continuum limit n —» oo with fixing // = 0(1). This means that the temperature approaches 
to the critical value as K — K* = 0(n -1 ) in taking the continuum limit. 

In this parametrization, 7 / can be expressed as 

C08h7 ' = 2 ^ + 1+2 ™ 2 ^' (3 ' 2) 
Since we are interested in the continuum limit, only the region / < n is relevant. Then, in 
such a region, 77 can be expanded for n> 1 as 

= +oin " ) ' (3 ' 3) 
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Substituting it to Eq.(2.28) and taking the limit of n —>• oo, Pi, - ■ • , Pi become 

2 


lim P\ = T 

n —Voo ■*- 


k =0 


lim P -2 = T 

n—>oo 

fc =0 

lim -P 3 = T 

n—>-00 - 1 - 

fc=l 

00 

lim P 4 = T" 

n — 


27 r *(fc+i) (r 1 +ir 2 yi+(^^) 2 ) 
1 + e V / 

^ 27ri(fc+i)^ri+ir 2 ^/ 1 +(( 5 ^ r ) 2 j 

1 + g 27rife ( Ti + iT2 \/i+(if;;) 2 ) 
1 _ e 27 rifc ( T i+ iT 2 V / l+(^*:) 2 ) 


(3.4) 


k =1 

where we have used the reflection property of 7 / (2.24). We can also show 

lim e~ mi ° = e _T2/i , 

n—>• 00 

and 


f 7/ = e 7 rT V(i) 2 + (^) 2 _ * 2 ri/ 


lim e 2 ,J = e 

n—> 00 


(3.5) 


(3.6) 


Combining Eqs.(3.4), (3.5) and (3.6), we see that the partition function in the continuum 
limit can be expressed as 


^cont. 
J i 1 


z = cJ 2 Zi 

i— 1 

where C is an irrelevant constant and Zf ont ' = lim n _ i . 0O Zj are given by 

jj7ri(fc+i) ^ri+ir 2 yi+(^^rj 


(3.7) 


z c° nt . = j || e irr 2 ^(fc+i ) 2 + (^) 2 j || 


^fcez 


fc =0 


1 + e 


Z cont. _ 
2 — 


/ 7 Cont. _ 
— 


Z cont. _ 
4 — 


e W(fc+i ) 2 + (£) 2N \ J] 


^fcez 


^fcez 


^fcez 


k=0 


^ 2 iri(fc+i) ^ri+mi+ (' j 


(3.8) 


OO 

(1 + e-^) n 

fe=i 

1 + e 27 r * fc ( Tl+ * T 2 V /l+ (^) 2 J 

00 

(1 - e-q n 

fe=l 

1 _ e 2 ™ fc ( T i+*W 1 + ( 3 ^) 2 ) 


which exactly agree* with the result in Eq.(BlO) of continuum field theory of a free massive 
Major ana fermion 

V'COnt. _ ry V’COnt. _ ry r^COnt. _ ry ryCOIit. _ ry (O C\\ 

Zi — Z 11 , z 2 — z 0j i, z 3 — Zi j0 , z 4 — zo, 0 ) 13.yj 


! The prefactor in E.(2.25) gives a non-universal overall factor in the continuum limit. 
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where the labels /x, u = 0, \ of Z^ v specify the boundary condition of the fermion as described 
in Appendix B. This means that the continuum limit of the two-dimensional Ising model 
with the boundary condition in Eq.(2.2) is the two-dimensional massive fermion theory on 
the torus. In particular, the combination ^ + iff becomes the complex structure r = T\ + IT 2 
of the continuum torus as shown in Eq.(2.29), and the parameter g in Eq.(3.1) representing 
the deviation from the critical temperature is nothing but the mass parameter of the fermion 
in the continuum limit. 

Before closing this section, it is worth looking at the RG structure of the Ising model. Fixing 
parameters t\ and 72 in Eq.(2.29) which specify the geometry of the torus, the parameter 
space of the Ising model is spanned by (n, g), and the procedure of taking the continuum limit 
is nothing but taking the limit of n —l 00 with fixing g. As we have seen, after taking this 
limit, the parameter g is precisely identical to the mass parameter of the two-dimensional 
massive free fermion and the limit g 0 corresponds to the massless Majorana fermion, 
that is, c = | CFT. Therefore the flow parametrized by the parameter g starting from the 
conformal fixed point exactly equals to the mass deformation of the c = ^ CFT (Fig. 2). 



Fig. 2 The parameter space of the 2D Ising model with the condition p = T\ n and m = T2n. 
The curved lines express the RG flows of an irrelevant operator associated with A After 
taking the continuum limit, the parameter g is identical to the mass of the two-dimensional 
free fermion. 


4. Holographic description of the RG flow 

4-1. Holographic RG flow of single scalar field 

In this section, we holographically describe the RG flow starting from the conformal fixed 
point along the parameter g. 

In the spirit of the AdS/CFT correspondence, the source (coupling constant) of an operator 
in the boundary field theory is identified with a field in the bulk gravity. The value of the 
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coupling constant of the boundary field theory varies towards infrared (IR) as a RG flow. On 
the other hand, the classical solution of the bulk gravity provides a trajectory of the bulk 
field along the radial direction of an asymptotic AdS geometry. This trajectory is regarded 
as the RG flow of the coupling constant of the boundary field theory away from a CFT at 
ultraviolet (UV) fixed point. As mentioned in the introduction, the two-dimensional c = ^ 
CFT is conjectured to be dual to the “quantum” pure gravity under the assumption that the 
path integral over the metric of 3D space-time is localized to the classical solutions, that is, 
the BTZ black holes. We here assume that the same is true at least in the neighborhood of 
the CFT fixed point; the classical solutions are dominant even in the presence of additional 
fields in evaluating the partition function of the quantum gravity. 

In the analysis of the previous section, we found that there are two independent flows from 
the conformal fixed point parametrized by g and 1/n. We expect that the parameter 1/n 
in the boundary theory corresponds to certain discretized version of bulk gravity, which is 
difficult to work out at present. In the following, we explore the continuum bulk gravity to 
study the RG flow of the parameter g, which can be identified as the mass of the fermion of 
the boundary field theory. 

Since the parameter g couples to the operator TT at the boundary field theory, it is 
plausible to assume that the corresponding field in the gravity is a real scalar field 0, which 
may be considered as the minimum number of degrees of freedom to describe the RG flow of 
a single parameter g. We thus consider the following action of three dimensional Euclidean 
gravity with a real scalar field 0; t 


S = 


1 


J d 3 Xyftj 


-R + R(0) + ^iF(0)d M 0cU 


16 ttG N ' ~~ V " " ' ' ^ ’ t 4 ' 1 ) 

where Gjy is the Newton constant of three-dimensional gravity, (g,u = 1,2,3) is the 
metric, R is the Ricci scalar, g = det(g^j,), and IF(0) is a function of 0, which describes the 
nonlinearity of the kinetic term of scalar held 0. Since the AdS 3 geometry with the radius 
L should be a solution of this system when 0 = 0, we demand V (0) to satisfy 


= 0 ) = 


(4.2) 


Since 0 should have a nonsingular kinetic term at least for small 0, we further require that 
K(cj>) is regular at 0 = 0. Without loss of generality, we can fix 


K(tj> = 0) = 1 


(4.3) 


by choosing the normalization of the held 0. Our hnal task is to determine functions V (0) 
and AT(0) by requiring the solutions 0,5^,, to describe the holographic RG how of the Ising 
model off critical temperature. 


1 We here simply omit to write the boundary terms of the gravity action. 
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4-2. Solution of the gravity corresponding to the mass deformation 

Since we are interested in the evolution of the scalar field along the radial direction of an 

asymptotic AdS geometry, we set the following ansatz for the metric and the scalar field: 

2 

ds 1 2 = e 2h ^dr 2 + ^ dx 2 , <f> = 4>(r), (4.4) 

*= 1,2 

where Xi {i = 1,2) express the two-dimensional transverse directions, and r € [ 0 ,oo) is the 
radial coordinate, which may be regarded as the Euclidean time. We have also assumed that 
the functions h(r) and </>(r) depend only on r. Note that we have fixed the gauge by setting 
9ij = j 2 dij in Eq.(4.4). Since the geometry is asymptotically anti de Sitter space, e 2h must 
be expandable around r = oo as 

e 2h = ^ + 0 (r- 4 ). (4.5) 

In this setup, the independent field equations are given by-f 

l -K<mr) 2 = 4 + hd, (4.6) 

VW = (-1 + hh) , (4.7) 

where the dot ' expresses the derivation with respect to r. Our first task is to obtain h{r ) 
and <p(r) for given fixed functions K(<f) and V(f) by solving Eqs. (4.6) and (4.7). Using the 
boundary condition in Eq.(4.5), we can integrate Eq.(4.6) to obtain h(r) in terms of <f(r ) 
and K(cj)) as 

e 2 h (r) _ f±_ e -sK(cf>(s))(<j>(s)) 2 ds' 

The standard way to solve the equations would be plugging this to Eq.(4.7) and solve the 
obtained single equation of <f(r ) for given K((f>) and V(f)- However we proceed the opposite 
way: We first fix the behavior of <f(r) from a physical requirement and determine the relation 
between K(tf) and V(4>). 

Since the mass deformation keeps the theory free, the scaling of the mass parameter should 
be trivial. Then we can assume 

h( a ) = —ho, (4.9) 

a 

where a is the length scale of the (boundary) field theory and go is the mass at the ref¬ 
erence length scale ao■ We should also recall that, in the context of the holographic RG, 
the radial coordinate r is identified with the RG parameter, namely the “Euclidean time”. 
The evolution of a bulk field along this radial direction is identified with the RG flow of the 
corresponding coupling constant. In our choice of gauge in Eq.(4.4), the radial coordinate r 


1 The field equation of the scalar matter (f> can be derived from those of the metric. 
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can be identified as the length scale a of the two-dimensional boundary held theory. Then 
the behavior of the held <f(r) can be regarded as the scaling behavior of the corresponding 
coupling constant [37, 38]. In the present context, we are looking for such a solution that 
corresponds to the mass parameter fi of the free fermion theory which exactly behaves as 
Eq.(4.9). To achieve this goal, we should hx the solution of the scalar held as a function of 
radial coordinate r (length scale parameter) 

fir) = -</> 0 = cj)(r) (4.10) 

r 

with a constant fo, instead of solving fir) for given hxed K{f) and V(cj>). Eq.(4.10) allows 
us to change a variable from r to f, and Eq.(4.8) becomes 

g2 Hr) _ £_ e - fo (r) x R ( x )dx' (4.11) 

Substituting Eqs.(4.10) and (4.11) into Eq.(4.7), we obtain a relation between Kif) and 

V(<f>) 

vif) = -ji (i - efo * xKix)dx • ( 4 - 12 ) 

This is a necessary condition for the bulk scalar held f to become the holographic dual to 
the mass parameter of the free fermion. 


4-3. Holographic RG flow in terms of Hamilton-Jacobi equation 

Even after imposing the exact scaling behavior of the mass parameter in Eqs.(4.9) and (4.10), 
we still have a freedom to choose one function of f in the action, either K{f) or V(cj>), in 
order to represent the holographic RG how of our system. In fact, with the requirement of the 
scalar held behavior fir) in Eq.(4.10), the bulk geometry is uniquely determined by Eq.(4.8) 
for any choice of either Kif) or Vif) which are related by Eq.(4.12) as a consequence of 
held equations. 

In order to hx it, let us further consider the holographic RG structure of this system based 
on the Hamilton-Jacobi equation of the bulk gravity[37, 38]: 

2 


G. 


ij:kl 


( 1 ss cl \ 


+ 


1 


W$<W V y/9 bgki J 2K((f>) \y/g 5(f) 


1 SS d 


= V(f) -R + (4.13) 


where Gij-ki is dehned by Gij-ki = gikdji ~ 9ij9kh an d R is the two-dimensional scalar cur¬ 
vature constructed from gij(x). The classical action S'd = S c \[gij(x), f(x)\ is a functional of 
the boundary values of the metric gij(x) ( i,j = 1,2) and the scalar held 4>(x) on the two- 
dimensional surface at a specihc value of the radial coordinate, say r = ro, and is obtained 
by substituting the classical solution into the bulk action in Eq.(4.1). 

The momentum constraint of the bulk gravity insures that the classical action is invariant 
under the diffeomorphism of the two-dimensional boundary. We can then expand the classical 
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action S c \ in powers of derivatives^ in the transverse directions: 

S cl = f d 2 x^{W^) + •••), (4.14) 

where • • • includes terms with the derivatives of <f> and the two-dimensional curvature tensors, 
which vanish when <fi is independent of x 1 and the transverse geometry is flat. The Hamilton- 
Jacobi equation in Eq.(4.13) can also be expanded in powers of derivatives, and gives the 
following relation among W((f>), K{(j)) and V((f >) at the leading order of the expansion in 
powers of derivatives[37, 38] 

-W'(4,f. (4.15) 


n*)—k w (*r + w[*) 


By combining this relation with Eq.(4.12), we can determine W(4>) in terms of K(cf>) as 

xK(x)dx\ . (4-16) 


2 /1 n 

WM = I <xp[ 1 j o / 

Also, repeating the argument given in Refs. [37, 38], we can obtain the /3-function of <f> as 

2 W'(<l>) 


P{4>) = 


K(cf>) W(d>) ’ 


and the holographic c-function 


c( 0 ) = 


1 


(4-17) 


(4.18) 


G n W(cf>) ’ 

which is a monotonically decreasing function when the coefficient of the kinetic term K((f>) 
is positive definite. Note that the c-function (4.18) is defined as the coefficient of the scalar 
curvature appearing in evaluating the expectation value of the energy-momentum tensor. 
This kind of the c-function has also been proposed using functional renormalization group 
equation [48]. 

The gauge/gravity correspondence asserts that the classical action of the bulk gravity is 
regarded as the (regularized) free energy of the dual boundary field theory. In our case, we 
have considered the three-dimensional gravity with a single bulk scalar field (j) as a candi¬ 
date of the dual description of the massive two-dimensional free fermion. Therefore one may 
naively expect that e~ 5cl should be equal to the partition function of the two-dimensional 
massive free fermion in Eq.(3.7), when </> is independent of the transverse coordinates x l 
and the transverse geometry is flat. However, we should recall that the Ising held the¬ 
ory corresponds to the three-dimensional “quantum” gravity [5, 6 ]. The three-dimensional 
gravity with appropriate boundary conditions possesses the conformal symmetry with the 
central charge | at the boundary[4]. In Ref. [27], as mentioned in Introduction, the three- 
dimensional pure gravity is considered. It has been show that the integration over all the 


§ In Refs. [37, 38], the authors divide the classical action into the local and the non-local parts as 
S[<f>(x), gij{x)] = S\oc[4>{x),gij{x)] + T[cj)(x), gij(x)], and show that the non-local part T satisfies the 
RG equation of the boundary field theory. In this paper, we further expand T[</>(x),gij(x)\ in powers 
of derivatives. 
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possible three-dimensional metric with fixing boundary condition is localized to the classical 
three-dimensional geometries (BTZ black holes) and the partition function of the quantum 
gravity turns out to that of the Ising field theory by relying on the boundary conformal 
symmetry. Our analysis in this paper is based on the same assumption that this quasi-semi- 
classical approach still makes sense even after adding the mass term to the boundary held 
theory. Therefore, the partition function of the massive fermion as a sum of functions of 
the moduli parameter r in Eq.(3.7) should be obtained after summing up all the possible 
classical geometries. 

However it is actually hard to carry out this procedure explicitly, since we look at only one 
solution of the classical held equations and the conformal symmetry is broken by introducing 
nonvanishing values of the scalar held <f> to represent the mass term of the boundary fermion. 
However it is still possible to determine the scalar potential V(4>) of the bulk gravity as 
follows. We observe that the boundary condition of the fermion becomes irrelevant and that 
the partition function becomes r-independent when the geometry of the boundary becomes 
R 2 . In this case, the partition function becomes extremely simple and we can expect that 
every solution gives the same contribution to the partition function even if there are several 
solutions in the bulk with the same boundary condition. When the geometry is R 2 , the free 
energy of the massive fermion is given by 

J7 = J Wf log ( p2 + ^ : (419) 

where /x is the mass parameter. In order for the free energy in Eq.(4.19) to be well-defined, 
we need to regularize both UV and IR divergences. Although the result depends on details 
of the regularization, the regularized free energy in general takes the form 

•T"reg = a + fyx 2 + c^x 2 log /r, (4.20) 

where a, b and c are some constants. 

As mentioned above, we regard the scalar field cj) at the boundary as the mass parameter 
of the boundary fermion and we identify the regularized free energy with the classical action 
S cl in Eq.(4.14) at the boundary which is obtained from the bulk three-dimensional gravity. 
In addition, W(<fi) in Eq.(4.14) is proportional to the free energy because cj) does not have 
^-dependence. Therefore can be expressed as 

W{(j)) = A + B(j) 2 + C(j) 2 log</>, (4.21) 

where A, B and C are some constants. Using Eqs.(4.16) and (4.15), we can determine AT(</>) 
and V(cj)) in the following. 

Suppose C / 0, we observe that IU(</>) is dominated by the last term C(f> 2 log 4> at small 
values of 4>. This implies that the RG flow at 4> —>• 0 is discontinuous with the AdS solution 
in Eq.(4.16) when C/ 0. Since we identify the evolution of </> from r = 0 to r = oo along the 
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radial direction (the Euclidean ’’time evolution”) as the RG flow from the UV fixed point 
to IR, this behavior is not acceptable as the RG flow of the mass parameter. Therefore, we 
should adopt a regularization scheme in the holographic renormalization group by imposing 
the condition 

c = C = 0. (4.22) 


Note that the dimensional regularization realizes this condition for example. With this choice 
of W(cp), we obtain K(<p) and W((pi) from Eq.(4.16) and the normalization condition Eq.(4.3), 


m) 


1 + 


and 

w M = I ( 1 + ■ 

Incidentally, we obtain the /3-function (4.17) as 


(4.23) 


(4.24) 


P(4>) = cp, 


(4.25) 


and the c-function from (4.18) as 

c (^) = TXZ' ( 4 - 26 ) 

The beta-function (4.25) is consistent with the notion that the bulk scalar field cp corre¬ 
sponds to the mass parameter of the boundary field theory, and the c-function (4.26) is a 
monotonically decreasing function of cp € [0, oo) from ^ to 0 as expected. Substituting this 
result into Eq.(4.12), we obtain 

VW = ~J2 (i + x) • (4 ' 27) 

which is monotonically decreasing V —>■ —oo as cp —>• oo (see fig.3). This is also consistent 
with the expectation that mass perturbation deforms the c = 1/2 Ising CFT to flow to a 
system with less degrees of freedom as dictated by the c-theorem. Namely it is likely to flow 
to nothing (c = 0) in the IR. 


5. Summary and Discussion 

In this article, we have obtained the partition function of the Ising model on Euclidean 
two-dimensional lattice with the twisted boundary condition representing the torus with a 
complex structure r in a discretized version. We have taken an appropriate scaling limit to 
obtain a continuum limit for the torus with the complex structure r, retaining a deviation 
from the critical temperature. The resulting continuum partition function agrees with that 
of the mass-deformed Ising CFT, namely the continuum field theory of massive Major ana 
fermion on the torus with the complex structure r. We have also discussed the RG flow 
of the Ising model off critical temperature in terms of the three-dimensional AdS gravity 
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Fig. 3 The shape of the potential corresponding to the flow from the c = \ CFT to the 
infrared by the mass deformation. The point cj) = 0 corresponds to the AdS boundary with 
the radius L and (f> = oo corresponds to r = 0. 


together with a scalar field cj) representing the RG flow of the mass parameter. By using 
the Hamilton-Jacobi equation, we have found that a simple nonlinear kinetic function R(</) 
in Eq.(4.23) and a simple potential V(cf>) in Eq.(4.27) for the scalar field </> describes the 
expected RG flow from c = 1/2 CFT at the UV region towards c = 0 at the IR. 

Let us comment on the relation of our results to those in Ref. [47] where a similar setup 
has been used to look for solutions interpolating two AdS geometries corresponding to two 
different conformal field theories. They assumed the canonical kinetic term K(cj>) = 1 for the 
scalar field cj), which is achieved by a field redefinition from our case. They required that the 
metric component Gtt has a double zero at the horizon ro (IR), at which the metric becomes 
another AdS geometry in addition to an AdS geometry at r = oo (UV). Near r —>• oo at UV, 
their solution exhibits a chirally asymmetric Virasoro algebra, where the excitation spectra 
of the left and right Virasoro algebras, L$ and Lq are different. This peculiar behavior comes 
from their Ansatz with the non-vanishing metric component Gap, where t, p are coordinates 
other than r corresponding to our X\,X 2 in Eq.(4.4). We have chosen a different ansatz with 
Gtip = 0 in Eq.(4.4) to obtain chirally symmetric Virasoro algebras, as in the Ising model. 

We have used the Hamilton-Jacobi equation to describe the RG flow of partition functions 
of the Ising model on the torus. There are more detailed informations such as multi-point 
correlation functions or (equivalently) partition functions on higher genus Riemann surfaces. 
In principle, one should be able to describe these data, which are worth studying. However, 
it is likely that one needs to overcome the problems associated with the strong coupling, or 
quantum effects in gravity in order to describe these data quantitatively in the presence of 
matter fields. 
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Recent applications of AdS/CFT to condensed matter physics provide many interesting 
insights into possible phase structures of strongly coupled system. In particular, introducing 
periodicity into the system to mimic lattice structures seems to be a crucial feature of 
these applications[51, 52, 53, 54, 55]. Since many realistic models of statistical physics or 
condensed matter physics are built on some discretized system, it is desirable to obtain 
discretized models as the boundary of some discretized model for the bulk in discussing the 
AdS/CFT correspondence. We hope that the continuum space-time would emerge by taking 
the continuum limit of the discretized system in the bulk and the gravity would be realized 
as a sort of a cooperative phenomenon. This ambitious objective is of course still quite hard 
to achieve. How to obtain a discretized version of (quantum) gravity realizing the given 
discretized model on the boundary is the key question which is worth pursuing. With this 
philosophy in mind and taking our result quite optimistic, we may expect that the usual 
two-dimensional Ising model would be a candidate of the discrete boundary model related 
to a discrete version of quantum gravity in the above sense, and our result may be a first 
step toward this direction, since the deviation parameter p from the critical temperature 
is actually included in the original discretized Ising model. However our discussion still 
remained within the continuum theory. In this sense, it would be a great step if one finds a 
gravity description of the RG flow parametrized by - in Fig. 2. Such an attempt is interesting 
in its own right, and furthermore can provide a concrete starting point to incorporate lattice 
structures into the AdS/CFT correspondence, which may play a vital role in condensed 
matter physics applications. 
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A. Transformation of H± and Sq- 

Let us express the generators of Spin(2n ) in the fundamental representation as 

\Jp,v) — i $p,a&vp) ■ (1 Y fl <C V Y 2 Tl) (Al) 

\ /per 

Let us consider H± and E-t defined in (2.16) in the fundamental representation, which are 
matrices with the size of 2n and we write them as H± and Sq, respectively, in the following. 


18/27 



They are written as 


H± = 


x y 
t 


0 


y 

0 y 


x y 
■t 


x 


■■■ 0 
0 ••• 0 

y ■■■ o 


E ± = a± 


\=F y 0 x ) 

where x and y are 2 x 2 matrices, 

I cosh 2 a cosh 26 i sinh 2 d cosh 26 \ 


/ 0 0 0 

1 2 0 0 

0 1 2 0 

\0 0 ••• 


••• =Fl2^ 
... o 
... o 

1 2 0 ) 


(A2) 


x = 


—i sinh 2 d cosh 26 cosh 2 a cosh 26 I 


( — i sinh 2 a sinh 26 —i sinh 2 a sinh 26 \ 

2 ~ • I 

i cosh a sinh 26 — | sinh 2 d sinh 26 ) 

(A3) 


and the overall factor a± in the definition of S-t is given by 

—i (n : even) 


a + = 


(A4) 


1 {n : even) 

, a_ = < 

—i (n : odd) 1 (n : odd) 

which are necessary in order to reproduce ( 2 . 8 ) in the spin representation. 

Our goal in this appendix is explicitly transforming the matrices H± and S-t into this 
canonical form. Then we can easily estimate Tr±H±Y,± in the spin representation by using 
( 2 . 21 ). 

We first introduce matrices fl-t € SO(2n) whose (i-j) blocks are given by 


i 

*’ y/n \ n 


/n \ n 


(i,j = !,-■■ ,n) (A5) 


,%3 y/H V n J v ~ n ° 

respectively, where R(0) is dehned in (2.20). The similarity transformations of H± and S-i- 
by become 

(M x N 2n „ i\ 

M 3 A^ 2n -3 


niH+ n . = 


n 3 

\N X 

(r(-D 


M 2 n—3 


M 2n -iJ 


(A 6 ) 




*(-?) 




Hi) / 
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and 


f Mo 


nt = 


^2n-2 0 \ 


m 4 


A^ n-4 


N 4 


M 2 n-4 


n 2 

0 


M2n-2 0 

0 MqJ 


(H~¥) 


Si T E_S!_ = 


B(-%) 


A(f) 


*(¥) 


V 

where Mj and Nj (/ = !,•••, 2n) are matrices, 


iCj 0 


with 


/ Aj iBA 

Mi = I , N I = i 

\—iBi Aj) VO —iCi, 


Aj = cosh 2 a cosh 2b — cos ( — j sinh 2 a sinh 26, 

\n J 

Bj = sinh 2a cosh 26 — cos f — ] cosh 2a sinh 26, 

V n J 

Ci = sin ( —) sinh 26. 

V n J 


(A7) 


(AS) 


(A9) 


Note that / runs odd (even) numbers for the matrices with the index + (—). It is easy to 
see that M/ and Ni satisfy 


M 2 n -i = Mi, N 2n -i = — Ni. (A10) 

Since A/, Bi and Ci satisfy 

Aj — B j — Cj = 1, (All) 

we can uniquely determine the parameters 7 / > 0 , 0 / G [ 0 , and e r = ±1 by 

Ai = cosh 7 /, Bi = e/ sinh 7 / cos 6 */, (7/ = ± sinh 77 sin 9i, (A12) 

where the sign in the definition of Ci takes + for 1 < I < n and — for n + 1 < / < 2n. 
Note that the 7 / appearing in (A12) is the same one defined in ( 2 . 22 ). We also note that 


20/27 







N 0 = N n = 0 and Mq and M n are given by 


M 0 = i?(±i7o), M n = 


(A13) 


where the sign appearing in the expression of Mq takes + in the disordered phase (a > b ) 
and — in the ordered phase (a < b). 

We can rearrange the matrices (A6) and (A7) by permuting the elements properly. To this 
end, we introduce the matrices which generate the transportations, 


1 ) 


T = 
± h3 ~ 


\ 


and the cyclic rotations, 


0 2 


L 2 (i-i-l) 


0 2 


*■2 (n-j)J 


, (1 < * < j < n) 


C'i,j — Tj—i,jTj—2 t j~i ■ ■ • Tij+ 1 - 


(A14) 


(A15) 


Then by defining 


S+ = 


S _ = 


^2,71^4,72 * ' ’ Cn—2,n 

(n : even) 

C n +i C 2 .-lC 4 .-l ' ' ' Ct7,_3 .— 1 

2 

(n : odd) 

C^ n _-^C 2 } n—2C4 5 7 t,_2 ' ' * C'n — 4 , 71—2 

(n : even) 

C2 5 72 — lC4 ?ri _l • • • Cfi —3,71 — 1 

(n : odd) 


(A16) 


we obtain 


(Q + S + ) T H + (n + S + ) = 


(n+s+y s + ( n + s + ) = 


(J) X 2r -1 (n : even) 

r =1 

re — 1 
2 

X 2r -i © M n (n : odd) 
(n : even) 


< r =1 


r =1 x 

(j)^ 4 ( 2r _ 1 7r ) © (-1 2 ) ) (n : odd) 
\ r =1 


n 


(A17) 
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and 


(ns-) T H- (n.S-) = < 


(f]_5_) r S_ (S4_S_) = < 


r =1 


© ^4 
k r =1 


n 


A 2r © M n © Afo 

(n : even) 

A 2r © Mq 

< r=1 

(n : odd) 

^-7T^ © (-1 2 ) © 12^ 

(n : even) 

VT J © 1 2 

(n : odd) 


(A18) 


where A"/ and P 4 (0) are 4x4 matrices defined by 


A/ = f ^ Ar M , P 4 (/>) = P(-</>) ©P(</>). 

\ N f M r I 


(A19) 


The matrices A j and P 4 ( <f>) are simultaneously transformed into the canonical forms as 
Pj A/P/ = P(-ze/7/) © P(-ie/7/), P/ T P 4 (^)P/ = R{~4>) © P(</>), (A20) 

using the 4x4 matrices, 

7 " 0 


/ cos4 


Pr = 


cos ■ 

0 
0 

—ej sin % 


cos ej sin ^ 0 


0 e/ sin 

h 
2 




\-6/ 


—e/ sin cos 

0 0 


0 


cos % 


(A21) 


t / 


Thus we define the matrices, 


2 



✓ n — 2 

2 


© 

£ 

1 

(n : even) 


P2r © I 4 

(n : even 

r= 1 

n— 1 

5 

P_ = < 

r=l 

n- 1 


P2r-1 © 12 

(n : odd) 


P 2r © 1 2 

(n : odd) 

r—1 



>. r=l 





We finally consider the combinations, 

t± = n±s±p ± , 


(A22) 


(A23) 


which transform in the fundamental representation into the canonical forms, respec¬ 

tively: 
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n: even 


Tl{H™± p + )T + 

= 0(- R (“* me 2r-l72r-l 
r= 1 


2 r - 1 


n 


2 r — 1 \ 

pn) © i?(-zme 2 r-i72r-i H- P 7r)), 

n / 


(—i) p | ^^R(-ime 2r 7 2r — —Ptt) © -R(—ime 2r 7 2r + — p7r)^ 

r=l 


R(-im^ n + pn) ® i?(±im7 0 ) j, 

(A24) 


n: odd 


rJ(i7^S^)T + 


= -Z 


0" {©(*(-< 


r =1 


-zme 2r _i7 2r ._i 


2 r — 1 


n 


— 1 \ 

pn) ® R(-ime2r-il2r-i H-P^r)) 

n / 

® R(—im^ n + p 7 r) , 


T T(H™£p_)T_ 


^^i?(-zme 2r 7 2r — ^P7T) ® i?(— ime2r r j2r + ~P 7r )) © A!(±im7o), 


r=l 


(A25) 


where we have used (A13). We can easily see det T± = 1. These results motivate to introduce 
(2.23). Note that, when we consider the continuum limit, we should approach to the critical 
temperature from the ordered phase. Thus the sign appearing in R(zLim'yo) is chosen as +. 

B. Partition function of 2D massive fermion on the torus 

Let us consider 2-torus with periods wi,w 2 €C and free Majorana fermion with mass M on 
it: 

S = \j d 2 x^ T D^>, (Bl) 

t\T 


(B2) 


where T is a two-component spinor, and D is the Dirac matrix given by 

D _ (di + id 2 M \ 
y — M d\ — i&2) 

In the following, we evaluate the partition function, 

Z = f d^e~ s = Pf D, (B3) 
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where Pf D denotes the Pfaftian of the Dirac operator D. 

Let {h,k 2 } = {—*072/^4, za;i/v4} denote the dual vectors corresponding to {wi,^}, where 
A = Im(o;2Wi) is the area of the torus. The plane-wave on the torus is given by 


( ~ — „2iri(~nRe(k 1 z)+mRe(k 2 z)) 

i\*i 6) — e 


_ i(,mu) 1 +nui 2 )z+i(muj 1 +nuj 2 )z] 


(B4) 


where z = x\ + 1 x 2 and z = x 1 — ix 2 and m, n take integer or half integer, m € Z + fi and 
neZ + i/, (/r, v = 0 or 1/2), depending on the periodicity of the fermion. Recall that the 
eigenvalue of the Laplacian —{d\ + d 2 ) = —Add is 


x ( 27r ) 2 | , ,2 

A mn = ^2 l mw l + nw 2| 

= A7 2 |m + rn| 2 , 

where A f = i 2 ™ 1 ^ and 

r = 0 J 2 /UJ 1 = T\ + ir 2 . (n, r 2 € 
Then the partition function in Eq.(B3) is expressed as 

z= Y 
0 , 1/2 


with 




[ N\J\m + nr + fi + vt\ 2 + ( M/M ) 2 . 


m,nEZ 

Rescaling the mass parameter M as 




27t M 

t 2 aa 


and using the zeta function regularization, are estimated as as 

>E„ e z V- 2 +(^) 2 (1 - e -^) ^ ! _ e 2 ™( Tl+iT2 \/i + (-) 2 ) 

n=l 

7 ri( 27 r+l) ( n+ira J 
1 — e V 


^0,0 = 6^' 


Z 0)i = e , rT 2 ^^ V ( n + i ) 2+ (^) 2 JJ 


n =0 


71 — p 71 ^ 2 

zi o — e 


V n2+ (^) 2 (1 + e -^) J ] 


n=l 


1 + e 


27 rm( T 1 -MT 2 


V 1+ (^) ; 


^11= g '"” 7 ’ 2 L„ e z \!{ n +\) +(£■) TT 

2 ’ 2 J. J. 

n =0 


00 7 ri( 2 n+l) ( Ti+ir 2 

1 + e 


\/ 1+ ((^) ) 


(B5) 


(B6) 


(B7) 


(B8) 


(B9) 


(BIO) 
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We should note that the exponents of the overall factors in Eq.(BlO) logarithmically diverge; 


7TT2 


7TT2 



T2/T 2 

2vr 


C(1 ) + otf), 


Tin 2 
4 7T 


C(1 ) + o(A 


(Bll) 


Thus, precisely speaking, we implicitly regard that the logarithmic divergence £(1) is properly 
regularized in Eq.(BlO). 
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